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Abstract 



About eight minutes before a black hole expires it has a decreasing mass of 10^° g, an 
increasing temperature of 1 TeV, and an increasing luminosity of 7 x 10^^ erg/s. I show 
that such a black hole is surrounded by a quasi-stationary shell of matter undergoing 
radial hydrodynamic expansion. The inner radius of this shell is bounded by ten times 
the Schwarzschild radius of 1.6 x 10~^ fm and has a temperature about one-tenth that of 
the black hole. The outer radius, as defined by the photosphere, is about 1000 fm, has a 
local temperature of 100 keV, and is moving with a Lorentz gamma factor of 10^ Most 
of the emitted radiation is in photons with small amounts in gravitons and neutrinos. 
I calculate the instantaneous photon spectrum and then integrate it over the last eight 
minutes to obtain the energy distribution dN^/dE = Airmp/lbE'^ for E > several TeV . 
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Hawking radiation from black holes |I| is of fundamental interest because it relies on 
the application of relativistic quantum field theory in the presence of the strong field limit 
of gravity, a so far unique situation. It is also of great interest because of the temperatures 
involved. A black hole with mass M radiates thermally with a temperature 
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where mp = G^^/^ = 1.22 x 10^^ GeV is the Planck mass. (I use units in which h = 
c = = 1.) In order for the black hole to evaporate it must have a temperature greater 
than that of the present-day black-body radiation of the universe of 2.7 K = 2.3x10"^ 
eV. This implies that M must be less than 0.8% of the mass of the Earth, hence the 
black hole must have been formed primordially and not from a stellar collapse. The black 
hole temperature eventually goes to infinity as its mass goes to zero, although once Th 
becomes comparable to the Planck mass the semi-classical calculation breaks down and the 
regime of full quantum gravity is entered. Only in two other situations are such enormous 
temperatures achievable: In the early universe (T similarly asymptotically high) and in 
central collisions of heavy nuclei like gold or lead (T = 500 MeV is expected at the RHIC 
(Relativistic Heavy Ion Collider) just completed at Brookhaven National Laboratory and 
T = 1 GeV is expected at the LHC (Large Hadron Collider) at CERN to be completed 
in 2005). The spontaneously broken chiral symmetry of QCD gets restored in a phase 
transition/rapid crossover at a temperature around 160 MeV while the spontaneously 
broken gauge symmetry in the electroweak sector of the standard model gets restored 
in a phase transition/rapid crossover at a temperature around 100 GeV. The fact that 
temperatures of the latter order of magnitude will never be achieved in a terrestrial 
experiment motivates me here to study the fate of primordial black holes during the final 
minutes of their lives when their temperatures have risen to 100 GeV and above. The 
fact that primordial black holes have not yet been observed does not deter me in the 
least. 

When Th <C rUe (electron mass) only photons, gravitons, and neutrinos will be created 
with any significant probability. These particles will not interact with each other but will 
be emitted into the surrounding space with the speed of light. Even when Th ~ me 
the Thomson cross section is too small to allow the photons to scatter very frequently 
in the rarified electron-positron plasma around the black hole. This may change when 
Th ~ 80—100 MeV when muons and charged pions are created in abundance. At somewhat 
higher temperatures hadrons are copiously produced and local thermal equilibrium may 
be achieved, although exactly how is an unsettled issue. Are hadrons emitted directly 
by the black hole? If so, they will be quite abundant at temperatures of order 150 MeV 
because their mass spectrum rises exponentially (Hagedorn growth as seen in the Particle 
Data Tables 0). Because they are so massive they move nonrelativistically and may 
form a very dense equilibrated gas around the black hole. But hadrons are composites of 
quarks and gluons, so perhaps quarks and gluon jets are emitted instead? These jets must 
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decay into the observable hadrons on a typical length scale of 1 fm and a typical time 
scale of 1 fm/c. Once the hadrons appear they may form an equilibrated gas around the 
black hole just as if they had been produced directly albeit with some time delay. One 
can find arguments both for and against ^ thermal equilibrium being maintained by 
the strongly interacting hadrons outside the Scharzschild radius. Certainly this is a very 
difficult and open problem in quantum statistical mechanics. Fortunately the answer is 
unambiguous for black hole temperatures greater than the electroweak scale of 100 GeV, 
as we shall see. 

I assume that a primordial black hole is surrounded by a shell of hydrodynamically 
expanding matter in local thermal equilibrium when Th > 100 GeV. This assumption will 
be shown to be self-consistent. A detailed description of how this situation comes to be 
is a difficult problem as discussed above and is not addressed in this paper. 

The relativistic hydrodynamic equations describing an adiabatic, steady-state, spher- 
ically symmetric flow with no net baryon number or electric charge and neglecting gravity 
were derived by Flammang and Thorne [Q. The equations may be cast in the following 
form. 

47rr^t>7^Ts(T) = Lin = luminosity (2) 
7T = constant (3) 

Here r is the radial coordinate, v is the local flow velocity with the associated Lorentz 
factor 7, T is the local temperature, and s(T) is the entropy density of the fluid. The 
Lin is the luminosity flowing into the fluid. It is generally less than the total luminosity 
of the black hole. For example, gravitons will escape without scattering from the matter 
and so will not participate in the hydrodyamic flow. These equations were derived in the 
context of accretion but, of course, they apply equally well to expansion. They were used 
to describe gamma-ray bursts from neutron stars by Paczyhski 0. Given the luminosity 
these equations allow one to flnd the functions v{r) and T(r) algebraically without solving 
any differential equations. 



A black hole has a Schwarzschild radius = 2M/m| = l/AnT. Note that 7rTh-2i? 



h 



1/2. Roughly, the average thermal momentum of a massless particle times the diameter 
of the black hole is 1/2. This is just a manifestation of the uncertainty principle applied 
to the creation of an excitation in a conflned region of space. The luminosity is 



dt = QAn^TMt (4) 

where a{M) is a function reflecting the species of particles available for creation in the 
gravitational fleld of the black hole. It is generally sufficient to consider only those parti- 
cles with mass less than Th; more massive particles are exponentially suppressed by the 
Boltzmann factor. I use 



a 



= 2.011 X 10"^ [4200Aro + 2035A^i/2 + 835Ni + gSiVa] . (5) 
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Here Ns is the net number of polarization degrees of freedom for all particles with spin 
s. The coefficients for spin 1/2, 1 and 2 were computed by Page and for spin by 
Sanchez 0. In the standard model A^o = 4 (Higgs), N1/2 = 90 (three generations of quarks 
and leptons), A^i = 24 (SU(3)xSU(2)xU(l) gauge theory), and N2 = 2 (gravitons). 
This assumes Th is greater than the temperature for the electroweak gauge symmetry 
restoration [0. Numerically a{T^ > 100 GeV) = 4.43 x lO'^. 
The entropy of a black hole is given by the area formula. 

S, = mlnRl = 4vr^ = ^ (6) 

The entropy per unit time lost by the black hole —dS^/dt is to be equated with that flowing 
through the matter, 47rr^f 7s. Taking the ratio of this entropy rate to the luminosity allows 
us to determine the constant in the hydrodynamic equation. 

jT = Th (7) 

This result is independent of the function a(M). 

The entropy density of weakly interacting massless particles is 

siT) = ^d{T)T' (8) 

where d{T) is the number of bosonic degrees of freedom; fermions get counted with a 
weight of 7/8. Using this in eq. (2) it is found that f is a monotonically increasing 
function of r and that its minimum possible value is f min = 1/^3. This is a sensible value 
since that is the speed of sound of a weakly interacting gas of massless particles. The 
corresponding radius and fluid temperature are rmin = 1.43-Rh and T{r^i^) = 0.816Th. 
The former assumes that gravitons and neutrinos do not contribute to the entropy of 
the fluid. The actual radius at which the fluid may be considered as thermalized could 
very well be greater than rmin, but the fact that the latter is already greater than the 
Schwarzschild radius is eminently sensible. 

The expanding fluid is in local thermal equilibrium when the mean free paths or 
thermalization lengths of the particles are less than the characteristic distance over which 
the temperature changes. This characteristic distance is I = 'yT\dr /dT\ (the gamma takes 
into account the Lorentz transformation from the rest frame of the black hole to the rest 
frame of the fluid). When r is greater than a few times the Schwarzschild radius v is 
already approaching the speed of light and 7 1. Then r^T^ = 360a{T^) / ndiT) and 
neglecting the slow variation of d with T (except near a flrst order phase transition) 

/ - ^ 
K = Q 
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Above 100 GeV the numerical value of K (leaving out gravitons and neutrinos) is 0.069. 

In a relativistic plasma it is not sufficient to consider only 2-body reactions. Of 
prime importance in achieving and maintaining local thermal equilibrium are multi-body 
processes such as 2 — > 3 and 3 — 2, etc. This has been well-known when calculating 



quark-gluon plasma formation and evolution in high energy heavy ion collisions [|TT| and 
has been emphasized by Heckler in the context of black hole evaporation Q]. This is a 
formidable task in the standard model with its 16 species of particles. Instead I require 
that the Debye screening length for each of the gauge groups in the standard model be 
less than /. The Debye screening length is the inverse of the Debye screening mass 
where n = 1, 2, 3 for the gauge groups U(l), SU(2), SU(3). Generically m° oc QnT where 
Qn is the gauge coupling constant and the coefficient of proportionality is essentially the 
square root of the number of charge carriers |12|. For example, for color SU(3) 



m^=g,^l + N,lQT (10) 

where iVf is the number of light quark flavors at the temperature T. The numerical values 
of the gauge couplings are: gi = 0.344, g2 = 0.637, and g^ = 1.18 (evaluated at the scale 
mz) 0. So within a factor of about 2 I have mP ^ T. This estimate makes sense when 
compared to two other measures. The average energy of a massless particle is 3T and so 
is about three times the thermal DeBroglie wavelength. On the other hand Carrington 
and Kapusta calculated the mean time between two-body collisions in the standard 
model for temperatures greater than the electroweak symmetry restoration temperature 
in the process of calculating the viscosity in the relaxation time approximation. Averaged 
over all particle species in the standard model one may infer from that paper an average 
time of 3.7/T. Taking into account multi-body reactions would decrease that by about a 
factor of two to four. 

The conclusion to be drawn from the above is that local thermal equilibrium should 
be achieved when 1/T < I which translates to the condition T < Th/15. Once thermal 
equilibrium is achieved it is not lost as r increases because / increases as 1/T^ whereas 
the Debye screening length increases more slowly as 1/T. The smallest radius at which 
local equilibrium is maintained is obtained from the condition on the temperature to be 
approximately 13Rh- This means that the thermalized matter surrounding the black hole 
is far enough away so that it doesn't affect the Hawking radiation whose calculation had 
assumed that the black hole radiated into a vacuum. 

An amazing observation in experiments on central collisions of lead nuclei at a lab- 
oratory beam energy of 160 GeV per nucleon at the SPS at CERN is that the outgoing 
hadrons, such as nucleons, pions, and kaons, have spectra that appear to be very thermal 
|T^. In those collisions, quantities like local temperature, pressure, energy density, and 



particle density vary significantly over a length scale of several fm. For black hole tem- 
peratures greater than 100 GeV one may easily calculate / in that region of the expanding 
matter where it is in the hadronic phase, where the local temperature is on the order of 
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100 to 160 MeV, and find that it is much greater than several fm. Therefore thermal 
equilibrium is maintainted at least down to that range of temperatures. 

Following Paczyhski's analysis of neutron star gamma-ray bursters let us assume 
that the matter maintains thermal equilibrium all the way down to a local temperature 
less than the electron mass. Gravitons and neutrinos have already escaped, leaving only 
photons and a rarified plasma of nonrelativistic electrons and positrons. The photosphere 



is defined as that radius where the optical depth is unity [jTSl. Interaction of the photons 
with the electrons and positrons proceeds through the Thomson cross section ctt- Defining 
n± as the total density of non-relativistic electrons and positrons one has the Thomson 
mean free path for photons at temperature T. 

1 ?> I ml o /me \ , , , 

At = = 2 ^=34.1A -f exp(me/T) (11) 



Using d = 2 aX these temperatures, corresponding to the entropy being dominated by the 
photons, the condition At = I determines the temperature Tps of the photosphere. 



^|jexp(,VT,.) = ^ (12) 

Thus the photosphere temperature decreases as the black hole temperature increases, 
roughly in a logarithmic manner. Requiring Tpg < 100 keV means that Th > 680 GeV. 
For black hole temperatures greater than about 700 GeV essentially all of the luminosity 
is in the form of photons, with the gravitons emitted directly without rescattering, and 
the neutrinos having lost thermal contact at some very small radius and correspondingly 
high temperature. 

The Lorentz gamma factor at the photosphere is determined by eq. (7) to be 7ps = 
Th/Tpg. When the black hole temperature is 1 TeV the temperature of the photosphere 
is 100 keV and 7ps = 10^. Let us now calculate the radiation from a surface moving at 
relativistic velocity. The phase space density / is invariant. Hence 

/(p, X, t) = /(p', x', t') = ^ (13) 

exp(E7Tps) - 1 

where the unprimed variables refer to the black hole rest frame and the primed variables 
refer to the rest frame of the fluid. We desire the flux in the black hole rest frame. Using 
E' = 7ps-E(l — fps cos^) the spectral flux is 



dE 27r2 



d(cos 9) cos 9f{E, cos 9) . (14) 
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The denominator of / may be expanded in a Taylor series and integrated term by term. 
The resulting series cannot be summed in terms of elementary functions. However, in the 
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limit that 7ps ^ 1 and E ^ Tps/7ps (a very small energy indeed!) one gets 



h 



ln[l-exp(-E/2Th)] • (15) 



dE 27T^^l^ 

Integration over all energy gives the total flux (energy/area/time). 



One may also integrate eq. (13) exactly to get the total flux; in the limit 7ps ^ 1 one 
obtains the same result. The corresponding expressions for each flavor of neutrino are 



dF, E^Ty, 



dE 27r^r, 



ln[l + exp(-£;/2rh)] , 



7 2 7 2 

F = — -Y^T^ = — ^ (17) 

The temperature and Lorentz gamma-factor at which the neutrinos lose thermal contact 
will be addressed in another paper. 

Equating the total luminosity of the black hole, excluding that carried by gravitons 
and neutrinos, with 47rit!psF^ determines the radius of the photosphere. 

2 _ ISOaeff 

" ps 

Here aes = with a evaluated at > 100 GeV: the standard model degrees of 

freedom are assumed to hold all the way up to the Planck temperature. Starting with a 
black hole of temperature Th > 100 GeV the time it takes to evaporate/explode is 

At ^ ^ .3 . 19 

This is also the characteristic time scale for the rate of change of the luminosity of a 
black hole with temperature Th. It is many orders of magnitude greater than the time it 
takes a light signal to travel from the Schwarzschild radius out to the photosphere, thus 
justifying the quasi-stationary assumption. A black hole with temperature 1 TeV has a 
Schwarzschild radius of 1.57 x 10~^ fm, its photosphere is at Rpg — 970 fm, and has 464 
seconds to live. 

Knowing how the temperature, flow velocity and radius of the photosphere evolve 
with time allows us to determine the energy distribution of photons emitted during the 
final minutes of the black hole. Integrating from t — —At to t — gives 

dN^ 45 aesml f^/^'^i^ 



dE Att^ a E^ 



/ dxxMn(l-e-^). (20) 
Jo 
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Both this integrated distribution and the instantaneous one, AnR'^^dF^/dE, are insensitive 
to the precise location of the photosphere. In the hmit E » 2Th eq. (20) becomes 

dN^ 4:71 nip . . 

'dE^lE^' ^ ' 

This energy distribution should be valid for any primordial black hole for photon energies 
much larger than 1 TcV = 10^^ eV. 

The hot shell of matter surrounding a primordial black hole provides a theoretical test- 
ing ground rivaled only by the big bang itself. In addition to the questions already raised, 
one may contemplate baryon number violation at high temperature and how physics be- 
yond the standard model might be important in the last few minutes in the life of a 
primordial black hole. Experimental discovery of exploding black holes will be one of the 
great challenges in the new millenium. 
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